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We present forecast results for constraining the primordial non-Gaussianity from photometric 
surveys through a large-scale enhancement of the galaxy clustering amplitude. In photometric sur- 
veys, the distribution of observed galaxies at high redshifts suffers from the gravitational-lensing 
magnification, which systematically alters the number density for magnitude-limited galaxy sam- 
ples. We estimate size of the systematic bias in the best-fit cosmological parameters caused by the 
magnification effect, particularly focusing on the primordial non-Gaussianity. For upcoming deep 
and/or wide photometric surveys like HSC, DES and LSST, the best-fit value of the non-Gaussian 
parameter, /nl, obtained from the galaxy count data is highly biased, and the true values of /nl 
would typically go outside the 3-cr error of the biased confidence region, if we ignore the magni- 
fication effect in the theoretical template of angular power spectrum. The additional information 
from cosmic shear data helps not only to improve the constraint, but also to reduce the systematic 
bias. As a result, the size of systematic bias on /nl would become small enough compared to the 
expected 1-cr error for HSC and DES, but it would be still serious for deep surveys with z m > 1.5, 
like LSST. Tomographic technique improves the constraint on /nl by a factor of 2-3 compared to 
the one without tomography, but the systematic bias would increase. 
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I. INTRODUCTION 

Any hints on primordial non-Gaussianity would be 
fruitful to clarify the generation mechanism for primor- 
dial density fluctuations in the early stage of the Uni- 
verse. Since the single-field slow-roll inflationary scenario 
predicts nearly Gaussian fluctuations (e.g., [lHZl])' a de- 
tection of large primordial non-Gaussianity will rule out 
the simplest inflationary model and provides us a new 
insight into the physics in the early universe. 

Traditional and popular method to detect primordial 
non-Gaussianity is to measure the three-point correla- 
tions of statistical fields (e.g., This correlation 
vanishes in the Gaussian fields, and non- vanishing signals 
of the three-point correlations would provide information 
on primordial non-Gaussianity. 

On the other hand, recent numerical and theoretical 
studies (e.g., [Ill IHnliJ]) have revealed that the local- 
type non-Gaussianity, originating from the non-linear dy- 
namics of scalar fields on super horizon scales, can induce 
a large-scale enhancement in the galaxy clustering am- 
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plitude. In the local-type non-Gaussianity, the primor- 
dial fluctuations characterized by the Bardccn potential, 
<£(cc), are described by the Taylor expansion of Gaussian 
field </>(x) as (e.g., [IB III 13): 



$(z)=#z) + /nl(0 2 (z)-(0 2 )). 



(1) 



The non-vanishing parameter, /nl 0, implies a de- 
parture from the Gaussian statistics, and even a small 
value of /nl has been found to produce a scale-dependent 
galaxy bias, which is prominent on large scales and at 
high redshifts. With a help of this property, the con- 
straint on primordial non-Gaussianity has been obtained 
recently, combining photometric and spectroscopic sur- 
veys [2cH30( , and it turned out that the results are rather 
comparable to that from CMB observation. This con- 
straint mainly comes from the quasar data obtained from 
photometric surveys, which are wider and deeper than 
spectroscopic surveys with a limited observation time. In 
this respect, wide and deep photometric surveys planned 
to start in the near future such as Subaru Hyper Suprime- 
Cam (HSC) survey [H, Dark Energy Survey (DES) 11, 
and Large Synaptic Survey Telescope (LSST) survey (33() 
would provide a more stringent constraint on primordial 
non-Gaussianity. 

In those photometric surveys, the observed galaxy dis- 
tribution at high redshifts often suffers from the magni- 
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fication effect due to the gravitational lcnsing, which ap- 
parently changes the number density of observed galax- 
ies [34H361]. As increasing rcdshift, since the ampli- 
tude of the density fluctuations becomes small and the 
observed angular separation between any pairs of two 
sources decreases, the galaxy auto-angular power spec- 
trum is shifted to small scales with the amplitude de- 
crease. On the other hand, the lensing contribution on 
the angular correlations become significant at higher red- 
shifts. As a result, the contribution of the magnification 
effect on the galaxy auto correlations is expected to be 
significant not only at high redshifts, but also on large 
angular scales. We thus naively expect that the mag- 
nification effect can mimic the scale-dependent galaxy 
bias, and ignoring magnification effect in the theoretical 
template for angular correlations would lead to a biased 
estimation of primordial non-Gaussianity. 

In this paper, we study the potential impacts of 
the magnification effect on constraining primordial non- 
Gaussianity from the upcoming photometric surveys. 
There are several studies forecasting constraints on pri- 
mordial non-Gaussianity through the scale-dependent 
galaxy bias (e.g., [29|, l37H39j ). Here, we pay a partic- 
ular attention to the magnification effect, and quanti- 
tatively evaluate the systematic biases arising from the 
incorrect treatment of the magnification effect in estimat- 
ing the non-Gaussianity parameter /nl, which has never 
been considered in the previous forecast studies. Fur- 
ther, we discuss the role of the cross correlation statis- 
tics between galaxy number density and other observ- 
ables such as cosmic shear. This has been also never in- 
vestigated in previous works, since the cross correlation 
signals are basically insensitive to the primordial non- 
Gaussianity compared to the galaxy auto correlations. 
However, we found that the cosmic shear-galaxy count 
cross-correlations have large signal-to-noise ratios, and 
help not only to improve the constraint on /nl, but also 
to reduce the systematic bias. 

This paper is organized as follows. In Sec [Hi we briefly 
review the scale-dependent galaxy bias induced by the 
primordial non-Gaussianity, and describe how the mag- 
nification effect changes the number density of observed 
galaxies. We then give the formalism for the angular 
power spectra obtained from photometric galaxy surveys. 
In Sec lIIIl we explicitly compute the angular power spec- 
tra and calculate the signal-to-noise ratios for auto and 
cross power spectra of galaxy count and cosmic shear. In 
Sec lIVl based on the Fisher matrix formalism, we quan- 
titatively estimate the impact of magnification effect on 
the detection of primordial non-Gaussianity, particularly 
focusing on three representative surveys, i.e., HSC, DES 
and LSST. Finally, Scc|V]is devoted to the summary and 
discussion. 

Throughout the paper, all the angular power spectra 
are computed from the modified version of cosmological 
Boltzmann code, CAMB [4(|, with the following set of 
cosmological parameters assuming a flat Lambda-CDM 
model, which is consistent with WMAP7 results 41j ; 



n h h 2 = 0.022, fl m h 2 = 0.13, Q A = 0.72, n s = 0.96, 
A s = 2.4 x 10~ 9 , t = 0.086, and w = -1, for the den- 
sity parameters of baryon and matter, dark energy den- 
sity, scalar spectral index, scalar amplitude at k = 0.002 
Mpc -1 , reionization optical depth, dark-energy equation- 
of-state parameter, respectively. Unless otherwise stated, 
non-Gaussian parameter is set to /nl = 0. The non- 
linear power spectrum is computed according to the fit- 
ting formula given in Ref . [42j . 



II. PROBING PRIMORDIAL 
NON-GAUSSIANITY FROM PHOTOMETRIC 
SURVEYS 

A. Primordial non-Gaussianity imprinted on 
galaxy bias 

In the presence of local-type primordial non- 
Gaussianity, recent numerical and theoretical studies on 
the clustering of halos/galaxies (e.g., [Til - t^ ) show that 
there appears a scale-dependent enhancement of the clus- 
tering amplitude on very large scales. Theoretically, the 
scale-dependent property of the halo/galaxy bias can be 
explained by a tight correlation between long- wavelength 
and short- wavelength modes, which usually vanishes in 
the Gaussian case. Especially, the modulation of short- 
wavelength modes responsible for forming halos is in- 
duced by the Newton potential or Bardeen potential 
<&(x). This fact leads to a strong scale-dependence for 
the fluctuations of halo/galaxy number density on large 
scales, and in Fourier space, we obtain [2^, [43| 



g{k,z) = [b G + Ab{k,z)]8(k,z), 



(2) 



where g(k, z) and S(k, z) are the fluctuations of galaxy 
number density and matter density fluctuations, respec- 
tively, and the quantity 6q implies the galaxy bias in the 
case of Gaussian initial condition. The function A6(fc, z) 
represents the non-Gaussian correction, which is given by 
il 



Ab(k,z) = f Nh A 



3fl 



NG 



k 2 T{k)D(z) 



(3) 



Here, the quantity f2 m is the matter energy density, Ho 
denotes the Hubble parameter at present, D(z) is the 
linear growth rate, and T(k) is the transfer function for 
linear matter density fluctuations, which is set to unity in 
the limit k — > 0. Thus, in the large-scale limit (k — > 0), 
the second term in Eq.([3]) becomes dominant, and the 
enhancement of clustering amplitude is prominent in a 
scale-dependent way. Since this term is inversely propor- 
tional to the growth rate D(z), non-Gaussian correction 
becomes also significant at higher redshifts. 

Assuming the universality of mass function, the quan- 
tity, Ang, gives <5 c (6q — 1) 0, HI], and the quantity 
5 C = 1.68 is the critical density for a spherical collapse. 
As advocated by several papers, however, the quantity 
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Ang would not be simply related to the halo mass func- 
tion, but depends on the merger history of halo/galaxy 
samples p3^ . In other words, we may have to deter- 
mine Ang from the observational data in practice. If 
this is the case, the non-Gaussian parameter /nl would 
be completely degenerated with the quantity ^4ng, and 
we need additional information on /nl like the galaxy 
bispectrum in order to break the degeneracy. Our pri- 
mary focus here is to explore the impact of magnification 
effect on the non-Gaussian parameter, and we simply as- 
sume A NL = 5 c (b G — 1) in the subsequent analysis. The 
influence of the magnification effect is a generic issue to 
constrain /nl from the photometric surveys, and we ex- 
pect that the results in the paper are also applicable to 
the case to combine other observations in breaking the 
degeneracy with Ang- 

B. Magnification effect on galaxy number density 

The number density of galaxies obtained from photo- 
metric surveys often suffers from the magnification ef- 
fect by the weak gravitational lensing of the large-scale 
structure (e.g., |44T - |46j ). Since the gravitational lens- 
ing changes the apparent magnitude and the area of the 
patch in the observed sky, it also changes the observed 
galaxy number density. Denoting the zero-mean fluctu- 
ations of the observed galaxy number density along a 
direction 6 at redshift z by n(6, z), we have [5i| 

n(0, z) = g{6, z) + (5s(z) - 2)k(0, z), (4) 

where the quantity k(9, z) is the lensing convergence 
at the position of source galaxy, and characterizes the 
change of the size of images. The convergence is given 
by m 

k{9,z) = / d x p-r S(x9,X), (5) 

2 Jo Xi z ) 

with the function x( z ) being the comoving distance. The 
second (5s(z)k) and third term {—2k) in Eq. arise 
from the modification of apparent magnitude and the 
area of the patch in the observed sky by lensing, respec- 
tively. 

The magnitude of the lensing effect depends on the 
slope parameter, s(z). Denoting the number of galaxies 
at redshift z, brighter than the ma gnit ude to by N(z, < 
to), the quantity s(z) is defined by [46| 

s dlog 10 N(z,<m)^ 
dm 

Note that in addition to the correction in Eq. (j4|), there 
exists another possible contribution related to the lensing 
effect, which has been addressed in Refs. jUSl). That is, 
the observed galaxies are selected according not only to 
the magnitude cut, but also to the size cut, and the latter 
also alters the galaxy number density. Nevertheless, the 



effect of size cut is basically proportional to the lensing 
convergence, and can be incorporated into the expression 
((4]), with a slight change of the meaning of slope index, 
s(z). In this respect, the results in the present paper is 
general, and applicable to the case taking account of the 
size cut. 



C. The angular power spectra 

The angular power spectra are the fundamental statis- 
tical quantity obtained from the photometric survey, and 
have a rich cosmological information. Here we write down 
the expressions for angular power spectra of galaxy num- 
ber counts and cosmic shear obtained from photometric 
surveys. 

The galaxy number density observed via photometric 
survey is projected onto the two-dimensional sky, and 
redshift information is obtained by dividing photometric 
galaxy samples into several subsamples binned with red- 
shifts. With the redshift distribution of galaxies in i-th 
bin, Ni(z), the two-dimensional distribution of observed 
galaxies in i-th bin, n,(0), are given by 

m{0)= f dz^-n(6,z). (7) 

The quantity Ni is the average number density per 
square arcminute in z-th bin, denned by 

poo 

Ni= / dz s Ni(z s ). (8) 
Jo 

On the other hand, the cosmic shear is measured from the 
cllipticity of each galaxy image. Using the photometric 
redshift information, we can also divide the estimated 
shear into several redshift bins. We denote the cosmic 
shear field in the i-th redshift bin by Ji(0). Then, the 
angular power spectra between the observables, X and Y 
(X and Y are either of 7.; or rij ) are given by the following 
expression: 

C f Y = lj% ^it(fc)Af (fc)Af (fc), (9) 

where Pi n it(k) is the matter power spectrum at an early 
time and k is the Fourier wave number. The functions 
Af (&) and Aj (k) are one of the following [H-[5l|]: 

AT(k) = ^±^V e (k;N l (z)), (10) 

A^(fc) = fc 2 | dz[b G + Ab(k,z)]^-D(zMk X (z)) 
+ (5si-2)£(£ + l)Ve(k;Ni{z)). (11) 
The function jt. is the spherical Bessel function and 
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Ve(k; Ni(z)) is defined by 



V t {k;Ni{z)) 



3 S^rn Hr\ 



• Nj(z s ) 
dz s — = — 

Ni 



(12) 

The quantity Si is the slope index in the i-th redshift bin. 
Although the slope index seems to have a strong redshift 
dependence [HoL [52j , we here assume the constant slope 
index within each redshift bin, and study the effect of 
time varying slope index. 

Note that for the photometric redshift determination, 
the uncertainty arising from the photometric redshift er- 
ror is crucial for the cosmological analysis [53| . To mimic 
this effect, we suppose that the photometric redshift esti- 
mates arc distributed as a Gaussian with rms fluctuation 
<j(z). Then the actual redshift distribution for i-th galaxy 
subsamples over the range, Zi-i < z < Zi, is related to 
the redshift distribution of galaxies, N(z), as [54| 

1 , . r / — z\ ( Zi — z 

(13) 

where the function crfc(x) is the complementary error 
function defined by 



2 f°° 

erfc(x) = — = / dzexp(— z 2 ). 
V 71 " Jx 



(14) 



III. MAGNIFICATION EFFECT ON 
ANGULAR POWER SPECTRA AND 
SIGNAL-TO-NOISE RATIO 

In this section, adopting a simple model of time evo- 
lution for bias parameter 6qi we compute the angular 
power spectra, and show how the magnification effect 
changes the amplitude of angular power spectra. Also, 
based on the fiducial setup of future photometric sur- 
veys, we estimate the signal-to-noise ratios for auto- and 
cross-power spectra of galaxy counts and cosmic shear. 

In Eq. |J2J), while the parameter 6g is assumed to be 
scale-independent on large-scales, it would manifest a 
strong time dependence. We characterize this by intro- 
ducing the following function (e.g., [55U571 ]): 



ba = bn 



D( z y 



(15) 



Here, we set bo = 1.5 and b z = 2.0 for fiducial values 
of the galaxy bias parameters. Since the photometric 
redshift information is available for most of the photo- 
metric surveys, we employ the tomographic technique, 
and divide all the galaxy samples into the three redshift 
subsamples. The redshift ranges for each bin are chosen 
as, z < 0.739, 0.739 < z < 1.14, and 1.14 < z, so that 
each redshift bin has the same number of galaxies. We 



assume the redshift distribution of galaxies N(z) as (e.g., 



N(z) = N e 



3z 2 



2(0.64z m ) : 



■ exp 



0.64z n 



3/2 



(16) 



where the quantities N g and z m are the total number 
of galaxies per square arcminute and mean redshift . As 
typical values of the upcoming deep surveys, we set N s 
to 35 arcmin -2 and z m = 1.0. As for the photo- z error, 
we adopt the simple scaling relation [5~4| : 



a(z) = 0.03(1 + *). 



(17) 



Finally, the influence of magnification effect on the an- 
gular power spectra depends on the slope parameter, for 
which wc set si = 0.5, S2 = 1.0, and S3 = 1.5, close to the 
recently estimated values from the observations [HO, HH . 



A. Magnification effect on angular power spectra 

In the presence of the magnification effect, the galaxy 
auto-power spectra, C" ,nj , and the shear-galaxy cross- 
power spectra, CJ a% , can be separately decomposed into 



several pieces: 



'-'II O -|" On , 



(18) 
(19) 



where the subscripts </, and //.,• respectively represent the 
contribution of the pure galaxy clustering and magnifica- 
tion, which arc identified with the first and second terms 
in Eq. pip . That is, the power spectra involving gi and 
Hi can be computed by simply neglecting the second and 
first terms in Eq. ([TT|) . respectively. 

In Fig. [TJ the galaxy auto and shear-galaxy cross angu- 
lar power spectra (C™ 3 ™ 3 and C^ 3 ™ 3 ) are shown. We plot 
the power spectra in the presence/absence of the mag- 
nification effect for the Gaussian initial condition (black 
solid/red dotted), and also show the power spectra in the 
absence of the magnification effect for the non-Gaussian 
case, /nl = ±10 (magenta solid). Note that the power 
spectra plotted here are independent of N g in Eq. (jTSJ) 
(see Eqs. ©-02)). 

For the galaxy auto power spectra C™ 3 " 3 , there ap- 
pears two contributions of the magnification effect in 
Eq. HID, i.e., C^ 3 and Cf 3 . As shown in the 
top panel of Fig. [TJ the cross power spectrum between 
galaxy counts and the magnification, C^ 3M3 (green dot- 
dashed), has similar angular dependence to Cf 333 (red 
dotted), and it slightly changes the overall amplitude of 
power spectra. The contribution of magnification auto- 
power spectra, C^ 3/i3 (blue dotted), has also similar £- 
dependence, but the amplitude exceeds Cf 3fl3 on large 
scales. This feature comes from the fact that the weak 
lensing effect is mainly attributed to the growth of struc- 
ture at lower rcdshifts. In the same manner, in the bot- 
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FIG. 1. The galaxy auto-power spectrum C™ n (top) and 
the shear-galaxy cross-power spectrum C] n (bottom). We 
plot the angular power spectra in the presence/absence of the 
magnification effect (black solid/red dotted) for the Gaussian 
initial condition, /nl = 0, and also show the angular power 
spectra in the absence of the magnification effect in the non- 
Gaussian case, /nl = ±10 (magenta solid). For comparison, 
we plot the contribution of the magnification effect (2*73/13, 
and /J3/13 in the top and 73/J3 in the bottom panel). 



torn panel of Fig. [TJ the amplitude of shear-galaxy cross- 
power spectra, Cj 393 , is enhanced at low-£ by the magni- 
fication effect. Hence, the magnification effect leads to a 
scale-dependent enhancement in the amplitude of angu- 
lar power spectra, which can mimic the effect of primor- 
dial non-Gaussianity through the scale-dependent galaxy 
bias. On the other hand, the bias parameters b and b z 
are independent of the scales, and their influences appear 
not only on large scales but also on small scales. Thus, 
even if treating these as free parameters and marginaliz- 
ing over them, the constraint on /nl seems to be unaf- 
fected by the galaxy bias parameters &o an d b z . 

In Fig. to elucidate the scale-dependent enhance- 
ment of the power spectra on large scales, we plot the 
fractional difference of the galaxy auto power spectra, 



FIG. 2. The dependence of galaxy auto-power spectra on the 
parameter /nl and the magnification effect for each tomo- 
graphic bin. We define the fractional difference of galaxy 
auto-power spectra, 5Ce/Ci which is given in two cases: 
Cf 9i (/NL = 10)/Cf 9i (/ NL = 0) - 1 (magenta lines) and 
C7'" 1 (/nl = 0)/Cf Si (/nl = 0) - 1 (black lines) . The solid, 
dashed and dotted lines represent = (1,1), (2,2) and 

(3, 3) bin, respectively. 



5Ce/Cg. Here, we examine the two cases: C 9i9i (f^h = 
10) /Cf 9i (f NL = 0) - 1 (magenta lines) and C™ 1 ™ 1 ('/nl = 
0)/C 9 ' 9 ' (/nl = 0) — 1 (black lines). As we expected, the 
impact of magnification effect is significant at higher red- 
shift bins. This is because, as increasing the source red- 
shifts, the gravitational lensing becomes significant and 
the amplitude of fluctuations g conversely decreases. The 
contribution of primordial non-Gaussianity is also signif- 
icant at higher redshifts, because the non-Gaussian cor- 
rection in the scale-dependent galaxy bias is proportional 
to the inverse of growth function, D(z) [see Eq.©]. Note, 
however, that the redshift dependence of the magnifica- 
tion effect (magenta lines) is rather different from that of 
the scale-dependent galaxy bias (black lines). This im- 
plies that the tomographic technique is useful to break 
the degeneracy between the effects of magnification and 
primordial non-Gaussianity. 



B. Signal-to- noise ratio 

Since the magnification effect enhances the amplitude 
of power spectra especially on large scales, the signal-to- 
noise ratio for angular power spectra would be changed. 
To estimate the size of this, in Fig. [3J we plot the signal- 
to-noise ratios for galaxy auto and shear-galaxy cross 
spectra in the presence or absence of the magnification 
effect. The signal-to-noise ratio, S/N, is defined by 
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IV. MAGNIFICATION EFFECT ON THE 
DETECTION OF PRIMORDIAL 
NON-GAUSSIANITY 

In this section, based on the Fisher matrix formalism, 
we now present the forecast constraint on /nl from pho- 
tometric surveys, and estimate the size of systematic bias 
from the incorrect treatment of magnification effect in the 
theoretical template of angular power spectra. As rep- 
resentative upcoming experiments for wide and/or deep 
surveys, we consider HSC experiment for wide but nar- 
row, DES for shallow but wide, and LSST for idealisti- 
cally deep and wide surveys. 



FIG. 3. The signal to noise ratio of each power spectra in 
the presence/absence of the magnification effect (red/green) 
for a survey with the mean redshift z m = 1.0 and the number 
density of galaxy N e = 35 arcmin -2 . We show the signal-to- 
noise ratio for all power spectra normalized by -J 1 / s k y since 
the signal-to-noise ratio is proportional to sj / sky . 



where the quantity AC^ y is the statistical error for each 
power spectra, given by 



A. Fisher matrix formalism 

Here, we summarize the Fisher matrix formalism used 
in the subsequent analysis, and describe the canonical 
survey setup for photometric galaxy surveys and CMB 
experiments. 

Given the angular power spectra theoretically 
parametrized by a set of parameters p, the Fisher ma- 
trix for the cosmological parameters is written as (e.g., 

M) 



(AC, 



XY} 



{21 + l)/ sky 



XY 



N, 



XY\ 



(21) 



Here, the parameter / s k y is the sky coverage of photo- 
metric survey, and Nf~ are the noise power spectra, 
which will be later given in Sec. IV [see Eqs. (|24l) - (|2"o]) ]. 
In Fig. [3J the signal-to-noise ratios are computed with 
Anax = 1000, and are normalized by / s k y . 

For the cross power spectra which have primarily no 
statistical correlation [i.e., C"'" 3 (i ^ j) and Cj^ 3 
(i < j)] [59j . the signal-to-noise ratios are significantly 
improved by the magnification effect. For example, the 
signal-to-noise ratio of CJ 2 ™ 3 increases by a factor of 
~ 115. This is because the magnification effect leads to 
a non- vanishing correlation between foreground galaxies 
and the background sources. On the other hand, the im- 
provement of the signal-to-noise ratio is relatively small 
for the power spectra which have strong statistical cor- 
relation even in the absence of the magnification effect 
[i.e., C"'™ J (i = j), and Cj 1 ™ 3 (i > j)]. Since the signals 
of these spectra are primarily very large enough to de- 
tect, the results indicate that the size of the constraint on 
/nl would not be drastically changed even if the magni- 
fication effect is properly taken into account in the data 
analysis. We finally note that, even in the absence of 
the magnification effect, the signal-to-noise ratios of the 
shear-galaxy cross correlations [C] tnj (i > j)] are com- 
parable to that of the galaxy auto correlations [C"'" J 
(i = j)]. This implies that, the parameter /nl would be 
constrained not only from the galaxy auto correlations, 
but also from the shear-galaxy cross correlations. 



F, 



y2l + l 

/ , o J sky 



1=1 



dpi 



dpj 



(22) 



where the quantity Ci represents the covariance matrix 
for the angular power spectra, pi is a cosmological pa- 
rameter which we want to estimate, and p fid is the set of 
fiducial cosmological parameters. 

In what follows, using a tomographic technique with 
the number of redshift bin, ./Vbin, we consider the num- 
ber density fluctuations ni, . . . , njv biI1 (or gi,...g^ hin 
when ignoring the magnification effect), and shear fields, 
71, . . . , 7jv bin , as observables obtained from the photo- 
metric surveys. Since these observables are rather sen- 
sitive to the late-time cosmic expansion and/or growth 
of structure, photometric surveys alone cannot give a 
tight constraint on all the cosmological parameters. To 
break the degeneracy between cosmological parameters 
and improve the constraints, wc include the informa- 
tion obtained from the primary CMB anisotropics by 
Planck [6l|. To be specific, we use the temperature (0) 
and (E-mode) polarization (E) data for primary CMB 
anisotropies. Denoting the noise power spectra by N^ Y , 
the full covariance matrix, Cp. is written as 



[Ci\ij - C 



XiXj 



N?' X3 5, 



11 ■■ 



(23) 



where Xi and Xj stand for m,...,n jv bta (or g\ , . . . , gN hin 
in the absence of the magnification effect), 71, ... , jN hin , 
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0, or E. The amplitude and shape of the noise spectra 
N* Y depends on the survey design which will be dis- 
cussed below. 

In the Fisher-matrix analysis, the forecast constraints 
depend on the properties of a photometric survey, charac- 
terized by the sky coverage, / s ky, the mean redshift, z m , 
and the total number of galaxies per square arcminute, 
N s . To show how the constraint on /nl depends on the 
survey design, we compute the Fisher matrix in the three 
representative photometric surveys; HSC for a deep sur- 
vey (/sky = 0.05, z m = 1.0 and iV g = 35 arcmin^ 2 ) and 
DES (/^y = 0.125, z m = 0.5 and N g = 12 arcmin- 2 ) for 
a wide imaging surveys, and LSST (/ s ky = 0.5, z m = 1.5 
and N g = 100 arcmin -2 ) as an idealistic survey, which 
is deeper and wider enough than the HSC and DES sur- 
veys. In Table HI we summarize the basic parameters of 
the survey design for three surveys used in the subsequent 
analysis. 

Let us now consider the noise spectra for each data set. 
In a photometric survey, the main noise source for galaxy 
counts is the shot noise given by 



N £ = Oij — , 

N 



(24) 



with the quantity N being the number density of galaxies 
per steradians in i-th redshift bin; 



Ni = 3600 N t 



(v 



str" 



(25) 



On the other hand, the noise source for cosmic shear 
measurement mainly comes from the intrinsic cllipticity 
of galaxies, which is described by 



N, 



1 /2 

The quantity (7 ; 2 nt ) is the rms intrinsic ellipticity. We 

adopt the empirically derived value, (Tint) 1 ^ 2 = 0.3 [Hf. 
In all surveys, the galaxy samples are divided into three 
redshift bins, and the ranges of redshift are chosen such 
that each redshift bin has same number of galaxies, N g /3. 
The resultant redshift ranges are summarized in Table HI1 
for each mean redshift z m . Note that the noise power 
spectra of CMB anisotropics, Nf e and Nf E , arc com- 



N 



(26) 



puted according to Eq. (3.3) of Ref.[63j, with the exper- 
imental specification for Planck [6lj . 

From the Fisher matrix, the 1-a (68 %C.L.) constraint 
on a cosmological parameter, cr(pi), marginalized over 
other parameters, is given by {F,^ 1 } 1 / 2 . Number of free 
parameters in the subsequent Fisher analysis is, in total, 
13, i.e., flbh 2 , ftmh 2 , Qa, n s , A s , r, w, /nl, in addition to 
the galaxy bias parameters (&o and b z ) and slope indices 
(s\„S2,Sz). Note that the number of free parameters is 
changed to 10 if we incorrectly neglect the magnification 
effect in the theoretical template of angular power spec- 
tra. The fiducial values of the parameters are the same 



TABLE I. Survey design for HSC, DES and LSST, i.e., the 
sky coverage, /sky, the mean redshift, z m , and the number of 
galaxies per square arcminute, Ng. In all surveys, the galaxy 
samples are divided into three redshift bin, and the ranges of 
redshift are chosen such that the each redshift bin has same 
number of galaxies, N g /3. The resultant redshift ranges are 
summarized in Table HIl for each mean redshift « m . 



Survey 


/sky Zm -*V g 


[arcmin a ] 


HSC [3LJ 


0.05 (2000deg ii ) 1.0 


35 


DES [32] 


0.125 (SOOOdeg^) 0.5 


12 


LSST [33 


0.5 (20000deg 2 ) 1.5 


100 



TABLE II. The relation between mean redshift, z m , and the 
redshift ranges of i-th bin computed in the case with Nbm = 3. 
Using Eq. (|13p . the redshift ranges are determined such that 
each redshift bin has same number of galaxies. 





redshift ranges 


0.5 


z < 0.369 0.369 < z < 0.569 0.569 < z 


1.0 


z < 0.739 0.739 < z < 1.14 1.14 < z 


1.5 


z < 1.11 1.11 < 2 < 1.72 1.72 < z 


2.0 


z < 1.48 1.48 < z < 2.29 2.29 < z 



as those in Sec lIIIi and we assume that the photo-z error 
and redshift distribution of galaxies are given in Eq. (|17p 
and Eq. pTJI) . respectively. 

Using the Fisher matrix formalism, we also evaluate 
the systematic bias in the best-fit value of cosmological 
parameters, Ap.;, arising from the incorrect treatment 
of the magnification effect in the theoretical template. 
Assuming the Gaussian likelihood function, the bias of 
the best-fit value, Api , can be estimated from 64 j 



^ 21+1 



sky 



x Tr lc- 1 (p)^-(p)C- 1 (p)(C(p) C(p)) 



p — pIlQ 

(27) 



where the covariance matrices C and C are computed 
with and without the magnification effect in the power 
spectra, respectively The Fisher matrix Fij is computed 
by ignoring the magnification in angular power spectra. 



B. Results 

1. Estimation of primordial non-Gaussianity for 
representative surveys 

Here we present forecast results for the constraint on 
primordial non-Gaussianity, especially focusing on the 
following cases: 

• +n: using galaxy counts alone taking into account 
the magnification effect in the theoretical template 

(cm 
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without magnifies I ion 
■■ +9 +9 + 1 

with magnification 




without magnification 
■■ +9 +9 + 1 

with magnification 
— +n +n + 7 
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2.4 1 , , r-i 2.4 
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-20 -10 10 20 




-20 -10 10 20 

2.4 n , , , ,-, 2.4 




-20 -10 10 20 

JNL 




1 1.5 

i>0 



FIG. 4. The 1-a error contours on &o-/nl (top left), & z -/nl 
(bottom left) and b z -bo (bottom right) planes for HSC. In each 
panel, we show the constraints in the case with (solid) and 
without (dashed) the magnification effect in theoretical tem- 
plate. In each case, we show the results obtained from galaxy 
number counts alone, and further including cosmic shear of 
galaxies. Note that CMB prior information from Planck is 
included in all cases. Information from the primary CMB 
is summed up to I = 3000, and other signals are included 
up to i'max = 1000. We also note that the fiducial values of 
the galaxy bias parameters and slope indices are chosen as 
b = 1.5, b z = 2.0, si = 0.5, s 2 = 1.0 and s 3 = 1.5. 



TABLE III. The 1-a constraints, £>"(/nl), and the system- 
atic bias, A/nl, on the primordial non-Gaussianity, in the 
case with/without the magnification effect in theoretical tem- 
plate. In each case, we show the resultant constraints ob- 
tained from galaxy number counts alone, and further includ- 
ing cosmic shear of galaxies. CMB prior information from 
Planck is included in all cases. We set the maximum multi- 
pole as £ma.x = 1000, but the primary CMB power spectra is 
used up to £ = 3000 in the Fisher matrix. 



FIG. 5. Same as Fig|4]but for DES. 





HSC 


DES 


LSST 




ct(/nl) 


A/nl 


ct(/nl) 


A/nl 


°"(/nl) 


A/nl 


+n 


4.8 




9.8 




0.86 




+9 


4.7 


11 


9.7 


21 


0.86 


7.1 


+n + 7 


3.5 




8.2 




0.49 




+5 + 7 


3.5 


0.19 


8.2 


1.7 


0.49 


-0.76 



• +g: using galaxy counts alone neglecting the mag- 
nification effect in the theoretical template (C 9 ' 9j ), 

• +n + 7: combining galaxy counts, cosmic shear, 
and their cross-correlations taking into account the 
effect of magnification (C™ ,nj , C^ n \ C^ l7j ), 

• +g + 7: combining galaxy counts, cosmic shear, 



1.6 

1.5 
1.4 


- (g) 


■J- 
® 


-2 2 

+; 



without magnification 
■■ +9 +9 + 7 

with magnification 
— +ti +n + 7 



-20 -10 10 20 



2.04 
2 

1.96 



2.04 
2 

1.96 






1.4 1.5 1.6 



-20 -10 10 20 

/nl 



1 1.5 

60 



FIG. 6. Same as Figg] but for LSST. The inset shows the 
zoomed error contours around the fiducial values. 



and their cross-correlations neglecting the effect of 
magnification, {Cf 9 \ Cf 93 , Cf l3 ). 

Note that, in all cases, we add the primary CMB infor- 
mation (i.e., Cf e , Cf E and Cf E ) in the Fisher matrix. 
The primary CMB power spectra are used to estimate 
the cosmological parameters up to I = 3000. In Table 
IIII1 we show the forecast results of 1-a constraints on 
/nl for three representative surveys. Also, in Figs|4]|6j 
two-dimensional contours of 1-a error on & z -/nl, ^o _ /nl 
and bo~b z planes are plotted in the cases of HSC, DES 
and LSST taking into account or neglecting the magnifi- 
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cation effect. 

Let us focus on the results from the galaxy counts 
alone, i.e., +n and +g (see red solid and dashed lines in 
FigsfJE]). Naively, the statistical error on /nl is expected 
to be large if wc properly take account of the magnifica- 
tion effect, because we need to specify the slope indices 
observationally and the number of parameters to be de- 
termined increases. As shown in FigsJUHl however, the 
statistical error on /nl does not change so much in all 
cases, and the fractional change is around a sub-percent 
level (see Table iHlj) . These figures also show that, as 
mentioned in Scc lIIII the degeneracy between /nl and 
the galaxy bias (bo and b z ) is weak (the correlation coef- 
ficient is ~ 0.1). On the other hand, the systematic bias 
on /nl arising from the incorrect treatment of the mag- 
nification effect is significant. As shown in Table. InTI the 
systematic bias is apparently very large for DES. Taking 
the ratio of the systematic bias to the statistical error, 
however, the largest value is obtained from the LSST 
case amongst three surveys. That is, the systematic bias 
is rather serious for LSST than for HSC or DES. 

Next discuss the importance of cosmic shear informa- 
tion, i.e., +n + 7 and +.9 + 7 (see green solid and dashed 
lines in Figs. |4]|6|). Compared to the cases with galaxy 
counts alone, the addition of cosmic shear data would 
not only improve the statistical error, but also reduce 
the systematic bias on /nl , irrespective of the treatment 
of the magnification effect. These results basically come 
from the non-vanishing shear-galaxy correlations, which 
also carry the information on /nl through the scale- 
dependent galaxy bias, as shown in Fig. [T] Table. IIIII 
shows that compared to the results with galaxy counts 
alone, the improvement of the constraint is by a factor 
of - 2 for LSST and of - 1.3 - 1.5 for HSC and DES. 
The size of systematic bias on /nl is now well within 
the 1-a statistical error for HST and DES, but it is still 
non-negligible for LSST. 



Dependence on maximum multipole ^ max , mean redshift 
z m , number of redshift bin Nhin, an d slope Si 
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FIG. 7. The systematic bias on /nl as a function of ^ ma x 
for +g (red) and +g + 7 (green) . The shaded regions denote 
the 1-a constraint on /nl for +n (thin) and +n + 7 (thick). 
Each panel show the case with z m = 0.5, 1.0, 1.5 and 2.0. 
The galaxy subsamples are divided into three redshift bins 
{Nhin. = 3). We assume N g = 100 arcmin -2 and / s k y = 0.5. 
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To elucidate the results in Sec II V B T1 in more details, 
we here consider the dependence of the systematic bias 
and statistical error on the various parameters, especially 
focusing on the LSST-like survey. Fig. [7J shows the sys- 
tematic bias A/nl using the tomographic technique with 
^Vbin = 3, plotted against the maximum multipole used 
in the parameter estimation, £ max . In each panel, the 
results obtained from a different survey depth are pre- 
sented; z m = 0.5, 1.0, 1.5 and 2.0 (in each case, red- 
shift ranges of each redshift bin are summarized in Ta- 
ble. HI]) . To infer the significance of the systematic bias, 
we also plot the 1-a constraint on /nl, c(/nl)i taking 
a proper account of the magnification effect (shaded re- 
gion). Fig. [5] also shows the same results as in Fig. [7J 
but, this time, we do not use the tomographic technique 
(i.e., -/Vhin = 1), and the slope index is simply set to 



FIG. 8. Same as Fig [7] but for TVbin = 1. In this case, we 
choose the slope as si = 1.0. Note that the range of y-axis is 
three times larger than that in Fig[7J 



Si = 1.0. Note that the systematic bias does not depend 
on /sky, while the 1-a constraint is proportional to U / s k y . 
Also, for a sufficiently large number density with N g s?10 
arcmin -2 , the results are almost insensitive to the choice 
of N s . Hence, the noise spectra are computed with fixed 
values of / s k y = 0.5 and N g = 100 arcmin -2 . 

Let us consider the results using the galaxy counts 
alone with iVbin = 3 (see FiglJj). As increasing z m , the 
statistical error on /nl (the thin shaded region) becomes 
rather improved, and the systematic bias (the red dashed 
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A(/nl) <t(/nl) 

Si 82 »s (Si, ss) = (0.5, 1.0, 1.5 + .' 

(0.5 + .s,1.0 ,1.5 ) +rj 

(0.5, 1.0 + a, 1.5 ) MM +n + 7 
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2„, = 1.5 
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tistical error on /nl- Fig. [9] shows the results in the case 
of the LSST survey. Here, the results are plotted against 
the shift of the slope index, s, defined by Si = s^fid + s - 
The lines with different color indicate the systematic bias 
of /nl obtained by varying the different slope index. Note 
that the statistical errors depicted as shaded regions are 
evaluated specifically in the case varying the slope S3. 

Fig. |9] indicates that the systematic bias of /nl is in- 
sensitive to the variation of slope index in lower redshift 
bin, but rather sensitive to it at higher redshift bin. This 
result is quite reasonable because the magnification ef- 
fect at high- z bin is significant compared to low-z bin 
(see Fig [2]). On the other hand, the statistical error on 
/nl is insensitive to the variation of slope index, as it is 
expected from Sec. IIII B[ In this respect, depending on 
the value of slope indices at high redshifts, the magnifi- 
cation effect on constraining primordial non-Gaussianity 
may become even more serious, and again, should be 
properly taken into account in the theoretical template. 



FIG. 9. Dependence of the systematic bias and 1-a constraint 
on the slope parameter. We vary only the one slope parameter 
(si, S2 or S3) and parametrized by a parameter s as sx = 
0.5 + s, S2 = 1.0 + s, or S3 = 1.5 + s. The parameter s is varied 
from —0.5 to 0.5. The red, green and blue lines show the 
systematic bias with varying ax, S2 and S3, respectively. The 
gray and black shaded regions represent the 1-a constraint 
in the case of +g and +g + 7, respectively, with varying the 
slope parameter in third redshift bin, S3. We assume LSST 
to calculate the noise spectrum. 



lines) is somehow reduced if we choose a large maximum 
multipolc. However, for a deep survey with z m Sri. 5 (bot- 
tom panels), the systematic bias is still non- negligible 
compared to the statistical error on /nl- This is true even 
if we include the cosmic shear data (the green dashed 
lines). 

Next discuss the importance of tomographic technique. 
Comparison between Figs. [5] and [7] implies that the to- 
mographic technique is quite helpful to reduce the sta- 
tistical error of /nl- Typically, with iVbin = 3, the error 
would be reduced by a factor of ~ 2 — 3. This is partly 
because the amplitude of the large-scale galaxy cluster- 
ing sensitively depends on the redshift in the presence of 
the scale-dependent galaxy bias (see Eq. (|3|) and Fig. [2J . 
Presumably, this may also break the degeneracy between 
/nl and slope indices. Note, however, that the tomo- 
graphic technique do not help to reduce the systematic 
bias, but rather it increases the size of bias. In this re- 
spect, a proper account of the magnification effect in the 
theoretical template is crucial in a LSST-like survey. 

So far, the fiducial values of the slope indices have been 
held fixed. In order to clarify the influence of the slope 
indices, keeping the number of redshift bins iVbin = 3, we 
allow to vary one of the three slope indices around the 
fiducial values, and estimate the systematic bias and sta- 



V. SUMMARY 

In this paper, we studied the impact of magnification 
effect on the detection of /nl from photometric survey. 
As representative upcoming photometric surveys, we con- 
sidered HSC for deep, and DES for wide, and LSST for 
an idealistically deep and wide survey. From the Fisher 
matrix analysis, we showed that, an incorrect treatment 
of the magnification effect on the theoretical template 
of angular power spectra leads to the systematic bias in 
the best-fit value of /nl ■ Especially, using galaxy counts 
alone, the size of systematic bias is significant for all three 
surveys (HSC, DES and LSST), and true values of /nl 
would typically go outside the 3-cr error of the biased 
confidence region. However, we found that additional in- 
formation from the cosmic shear observations helps not 
only to improve the constraint, but also to reduce the 
systematic bias on /nl- As a result, the systematic bias 
can become negligible for HSC and DES surveys, com- 
pared to their expected errors on /nl ■ A proper account 
of the magnification effect docs not increase the statis- 
tical error on /nl (~1% )■ Nevertheless, for LSST, a 
relative significance of the systematic bias still remains 
and the magnification effect should be correctly taken 
into account in the theoretical treatment. 

We further explored the various cases by changing 
parameters characterizing the survey properties, and 
showed that the tomographic technique using photomet- 
ric redshift information leads to a significant improve- 
ment on the statistical error on /nl, but it doe not help 
to reduce the systematic bias. In any case, high-z obser- 
vations arc indispensable for tightly constraining primor- 
dial non-Gaussianity, but the influence of the magnifica- 
tion effect would be inevitable. This is particularly true 
for deep imaging surveys like LSST (z m <rl.5). A proper 
account of the magnification effect in the theoretical tem- 
plate is thus quite essential for an unbiased estimate of 
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primordial non-Gaussianity. 
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